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ABSTRACT 


By using the Beam and Warming implicit-factored method of solution of the 
Nayier-Stokes equations, velocities were calculated inside axisymmetric piston- 
cylinder configurations during the intake and compression strokes. Results are 
presented in graphicial form which show the formation, growth and break-up of 
those vortices which form during the intake stroke by the jet issuing from the 
valve. It is shown that at bore-to-stroke ratio of less than unity, the vortices 
may break-up during the intake stroke. It is also shown that vortices which do 
not break-up during the intake stroke coalesce during the compression stroke. 
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I. INTRODUCTION 

Fluid motion inside cylinders of reciprocating piston engines has been 
studied both by experimental (1 - 12) and analyticial (13 - 34) methods. These 
investigations have revealed 1’ it recirculating flows (vortices) form during the 
intake stroke by the jet issuing from the intake valve. Under some conditions 
these vorticies may persist throughout the compression stroke (32) and seriously 
affect the ignition and flame propagation processes. Since these vortices can 
play such an important role in the combustion process* their formation, growth, 
and break-up during the intake stroke, and their subsequent behavior during the 
compression stroke must be understood. However, most of the previous studies 
focused only on the problem of vortex motion during the intake stroke. In spite 
of the importance of the problem, relatively little is known about the behavior 
of vortices during the compression stroke. This investigation was addressed 
thi-refore to the study of vortex motion during the compression stroke of those 
vortices which form during the intake stroke. To accomplish this, velocities 
were c’^leulated inside axisymmetric piston-cylinder configurations. The veloc- 
ity fields were presented in graphical forms which illustrate the formation, 
growth, and break-up of the vortices. 

II. DESCRIPTION OF THE PROBLEM 

The following problem was analyzed. A hollow circular cylinder is closed on 
one end by a flat piston and on the other end by a flat plate (Fig. 1). The 
piston is connected to a crank shaft through a connecting rod. The piston is 
driven by rotation of the crank shaft about the crank pin at an angular veloc- 
ity/1 resulting in a piston velocity u , 

r 
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The flat plate has a centrally located annular opening (valve opening) in it 
which opens instantaneously at the beginning of the intake stroke (crank angle t 
= 0) and closes instantaneously at the end of the intake stroke (^ « tr). 

The temperatures at the cylinder wall T^, valve T^, cylinder head Tj^, 
and piston T are constants, but may have different values. 

r 

The fluid enters the piston-cyl inder configuration described above through 
the valve opening during the intake stroke. The stagnation temperature T^- and 
stagnation pressure P. of the entering fluid are both taken to be constants. 

At the valve opening (seat angle a), the entering fluid may have velocity com- 
ponents in the radial (V^^) and axial (V^) directions, but not in the tangen- 
tial direction. The magnitude of the velocity at the valve opening depends upon 
the instantaneous flow field inside the cylinder as well as the stagnation tem- 
perature and pressure of the entering fluid. 

The viscous and thermally conducting fluid which enters the cylinder is an 
ideal gas having constant thermodynamic and transport properties. 

At the beginning of the intake stroke, tha gas in the clearance volume (re- 
sidual gas) is taken to be a^ stagnant, ideal gas at stagnation temperature T. 
and stagnation pressure P (P P.)» The residual gas has the same ther- 
modynamic and transport properties as the intake charge. 

III. METHOD OF SOLUTION 

The basic equations governing the problem are the conservation equations of 
mass, radial momentum, axial momentum, and energy (35). To render these equa- 
tions more amenable to numerical methods of solution, body forces were neglected 
and the bulk viscosity was taken to be zero. 
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The governing equations, applicable to laminar flows, are summarized in Ta- 
ble 1. Turbulence was taken into account by the method described subsequently. 
Equations 1 - 13 (Table 1) constitute a closed system in the four basic depen- 
dent variables: density (p), radial velocity (V^), axial velocity (V^), and 

energy (e). 

The boundary and initial conditions corresponding to Eqs. 1 - 13 are given 
in Table 2. The no-slip condition requires the fluid velocity next to a solid 
wall to be equal to the velocity of the solid wall (Eqs. 14 and 15). The gas 
temperature next to solid walls equals the temperature of the walls. With all 
walls maintained at constant temperatures, this boundary condition results in 
Eqs. 16 - 19. 

The stagnation pressure and stagnation temperature of the gas at the valve 
opening are constant with respect to time as reflected hy Eqs. 20 and 21. The 
radial and axial velocities at the valve opening are related, as indicated by 
Eq. 22. The axial velocity at the valve opening is determined by applying the 
conservation of mass equation at the valve opening (Eq. 23). 

The symmetry conditions at the center line result in Eq. 24. At .time t 
equal to zero, the piston is at the TDC position (^ = 0) and the residual gas in 
the clearance volume is a stagnant, ideal gas at stagnation pressure P and 
stagnation temperature T^. These initial conditions are *:pecified by Eq. 25. 

Recent evidence suggests that the type of turbulence model employed in the 
analysis does not affect significantly the calculated flow pattern (36). For 
this reason the increased mixing due to turbulence was simulated in a simple and 
convenient manner by choosing appropriate values of the effective transport 
Ifroperties. 
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It has been observed that the effective viscosity for turbulent flows inside 
spark ignition engine cylinders is roughly 100 times higher than the viscosity 
corresponding to lamitiar flows (37). Therefore, the viscosity, h, was taken to 
be 100 times the /,scos1ty of air. The valve of the thermal conductivity, a, 
was selected by taking the turbulent Prandtl number (P^ = MCp/x) to be equal 
to unity (38, 39), 

Solutions to the governing equations and the corresponding initial and 
boundary conditions formulated above must be obtained by numerical methods. In 
this investigation, the impl icit-factored, f inite-difr’erence method of Beam and 
Warming (40, 41) was used to obtain solutions. The details of the method are 
not given here because of the very long description this would require. Readers 
intereiited in the solution procedure are referred to references 42 and 43 which 
contain detailed descriptions of the method, including the formulation of the 
finite-difference equations and numerical boundary conditions, and discussions 
of the grid sizes, time steps, stability criteria, and computer time. 

IV. RESULTS 

Calculations were performed to explore the vortex patterns inside the cylin- 
der during the intake and compression strokes. The emphasis of this investiga- 
tion was on the behavior of vortices during the compression stroke. However, 
vortex formation during the intake stroke was also examined, since it is these 
vortices which carry over into the compression stroke. 

The ranges of the parameters for which numerical solutions were obtained are 
summarized in Table 3. The values of the parameters describing the geometry and 
operating conditions were selected so as to correspond to those typically found 
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in spark-ignition engines. The thermodynamic and transport properties were 
chosen so as to be physically reasonable ard to permit the use of convenient 
grid sizes in the numerical solutions. 

The axial and radial components of the gas velocities inside the cylinder 
were calculated as functions of time and for all the conditions listed in Table 
3. In the interest of brevity, only typical results are presented here which 
illustrate the major features of the flow pattern. The resultant velocity vec- 
tors were plotted providing a picture of the flow pattern. These graphical re- 
sults \<ere deemed adequate fov" the purpose of this study which was to observe 
the flow pattern inside the cylinder. 

Intake Stroke 

During the intake stroke, Vwo toroidal vortices were formed by the vorticity 
generated by the jet at the valve opening and by the adverse pressure gradients 
produced by the jet impinging on the piston surface (Figs, 2 and 3). One of 
these vortices was located between the jet and the cylinder wall (cylinder-head 
vortex). The other was located between the jet and the center line (valve vor- 
tex). The formation and subsequent motion of these vortices were found to be 
consistent with those observed in previous experimental (1 - 9) and numerical 
(15 - 30) studies. 

However, the present results showed one phenomenon that has not been re- 
ported previously. At bore-to-stroke ratio of less than unity, the cylinder- 
head vortex broke up into two smaller vortices (Fig. 2). The break-up was 
caused by interaction of the valve and cylinder-head vortices. Ekchian and 
Hoult (1) also noted vortex break-up in their tests with water injected into a 
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circular cylinder. However, in Ekchian and HouU's experiments, the vortex 
break-up was due to flow instability and not to vortex Interaction. In Ekchian 
and Hoult's experiments the vortex broke up and then degenerated Into a random 
flow. In the present study the vortex did not degenerate Into a random flow 
because flow instability and three-dimensional flow (needed for the complete 
vortex break-up) were not Included in the numerical solution. 

Compression Stroke 

During the intake stroke, the jet separated the cylinder-head and valve vor- 
tices minimizing the Interaction between them (Figs. 2 and 3). During the com- 
pression stroke, the cylinder-head and valve vortices were no longer separated 
by the jet, allowing the two vortices to Interact. Because of the rotational 
motion of the two vortices, they forced each other towards the piston. Since 
the two vortices were of unequal strength, the weaker one was pushed closer to 
the piston surface by the stronger one. 

As the piston speed increased during the compression stroke (crank angle ^ 
between 190 and 240 degrees), the two vortices began to coalesce (Fig. 4). The 
time of this coalescing depended upon the strength of the vortex which was 
nearer to the piston surface. 

At crank angles ^ between 230 and 240 degrees, the two vortices coalesced 
into a single toroidal vortex (Fig. 4). The coalescing of the cylinder-head and 
valve vortices has not been described by previous investigators. The disappear- 
ance of one vortex was found by Ashurst (18) and by Diwaker, (19) in nu- 

merical simulations of similar piston-cylinder problems. 
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Towards the end of the compression stroke (crank angle ^ > 320 degrees) a 
new recirculating flow (corner vortex) formed in the corner of the cylinder-head 
and cylinder wall (Fig. 4). The formation of the corner vortex has also been 
reported by Chong, et ^ (14) and Gosman, et (29). 

Finally, it is noted that in an actual cylinder, a vortex also forms near 
the cylinder wall due to the piston scraping off the boundary layer next to the 
cylinder wall during the compression stroke (10 - 12, 31, 32, 44 and 45). This 
vortex was not included in this study, since the interest here was only in those 
vortices which form during the intake stroke. 
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TABLE 1. GOVERNING EQUATIONS* 



L 


ORlGiNftt 
OF POOR 

TABLE 2. BOUNDARY AND INITIAL CONDITIONS 



Equation 


Eq. 1 

S'-O/'i-Mr 

piston surface 

14 


cylinder head, cylinder wall, and valve 

15 

T'‘Ty 

valve 

16 


cylinder head 

17 

rm%, 

cylinder wall 

18 


piston surface 

19 


valve opening 

20 


vavie opening 

21 

ot 

valve opening 

22 


valve opening 

23 

^my^m ■ 0 

center line 

24 

z«-R/(RTj/h),v^.-4-o 

RCCMC'r-R 

everywhere inside the piston-cylinder 
configuration at time t =* 0 

25 




TABLE 3, SUMMARY OF PARAMETERS USED IN THE NUMERICAL CALCULATIONS* 


Tp/rc 

(rpm) 

Pi/Pc 

a(degrees) 

0.78 

400 

1 

0 

1.67 

400 

1 

0 

1.67 

400 

1.036 

0 

1.67 

400 

1 

45 

1.67 

400 

1 

0 


The following parameters were the same for all cases studied: y « 3.2 x 
ig-3 kg/m-s, X * 3.224 kg-m/s2 - *K, h* = 42.252X/s2, Cp = 1,006.6 
m2/s2 - K, M » 28.96 kg/kg-mole, R = 8314.3 kg-m^/ kg-mole-*K-s, rp « 
0.05 m, ry « 0.01875 m, ry * 0.03438 m, 6(t=0) = 0.01 m, )q * 0.2m, Tj 
= Ty « Th » T^y » Tp = 340 K, and Pr = 101325 kg/m-s^. Results are 
shown only for the first two cases listed in this table. 
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